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Spontaneous pattern formation is a phenomenon ubiquitous in nature, examples ranging 
from Rayleigh-Benard convection to the emergence of complex organisms from a single cell. 
In physical systems, pattern formation is generally associated with the spontaneous break- 
ing of translation symmetry and is closely related to other symmetry-breaking phenomena, 
of which (anti-)ferromagnetism is a prominent example. Indeed, magnetic pattern forma- 
tion has been studied extensively in both solid-state materials and classical liquids. Here, 
we report on the spontaneous formation of wave-like magnetic patterns in a spinor Bose- 
Einstein condensate, extending those studies into the domain of quantum gases. We observe 
characteristic modes across a broad range of the magnetic field acting as a control parame- 
ter. Our measurements link pattern formation in these quantum systems to specific unstable 
modes obtainable from linear stability analysis. These investigations open new prospects for 
controlled studies of symmetry breaking and the appearance of structures in the quantum 
domain. 

Bose-Einstein condensates of alkali atoms such as ^^Rb offer unique opportunities to study 
classical non-linear effects in the quantum regime Governed by many-body quantum mechan- 
ics but subject only to weak interactions due to atomic collisions, Bose-Einstein condensates are 
very well described by a single-particle wavefunction whose dynamics is given by a non-linear 
Schrodinger equation known as Gross-Pitaevskii equation. Spinor Bose-Einstein condensates ^'^ 
in addition make use of the non-zero spin and associated magnetic moment of alkali atoms in 
their internal ground state. Adding the orientation of this spin as a degree of freedom, non-trivial 
dynamics arises from the spin-dependence of coUisional interactions on the one hand, and inter- 
action of the atomic moments with an external magnetic field on the other hand. Although based 
on fundamentally different interactions this mechanism effectively adds magnetic properties to the 
quantum gas, similar to the ones in solid state systems. 

Indeed the dominant spin dependent interaction term in ultracold atomic collisions is propor- 
tional to the product of hyperfine spins Fi • F2 and thus reminiscent of the Heisenberg model of 
solid-state magnetism - this analogy motivated the classification of spinor BEC ground states as 
ferromagnetic and anti-ferromagnetic in previous experiments However, it is worth noting 
that interactions of ultracold atoms are purely local (usually approximated by a 5-potential), in 
contrast to the next-neighbor interactions common in solid-state models. Further differences arise 
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from the hyperfine nature of the atomic spin, i.e. its composition of nuclear and electronic spins. 
First of all, its value can be significantly larger than the electron spin S' = 1/2, implying a higher 
dimensional state space with potentially much richer structure and dynamics. Secondly, the inter- 
nal structure of hyperfine states leads to a more complex interaction with external magnetic fields, 
beyond the linear Zeeman effect. The competition between spin-dependent coUisional interactions 
and quadratic Zeeman effect lies at the heart of spinor BEC physics. It is the driving force behind 
coherent oscillations in quasi-homogeneous systems and it is responsible for the quantum 
phase transition associated with the formation of irregular spatial spin structures 

The wavefunction of a spin-F condensate can be written in terms of a density and a local 
spinor, iprn{x) = v^n(x)Cm(^). normalized to the total number of particles / n{x) dx = N and 
unity J2m=-F iCmP = 1 respectively. In the simplest case of a spin-1 condensate, the local mean- 
field energy associated with is given by gou + gin{F)'^. The coefficients 5^0,1 are determined 
by the atomic s-wave scattering lengths. The gi term represents the spin-dependent part of the 
interactions, (F) is the local expectation value of the spin vector. This term generally is a small 
fraction of the total mean-field energy. For ^^Rb in the F = 2 state, gi is positive, indicating an 
anti-ferromagnetic ground state with (F) = 0. In contrast, for F = 1 gi is negative, favoring a 
maximally stretched state (ferromagnetic ground state). Further spin-dependent contributions arise 
from the interaction of the atomic momenta with an external magnetic field, — p(F^) + q{F^). The 
linear and quadratic Zeeman effect are parametrized hy p (x B and q oc B'^ respectively. In the 
absence of inhomogeneous fields, the linear Zeeman effect merely causes Larmor precession of 
the spin which disappears in a suitable co-rotating coordinate system. Consequently, all non-trivial 
dynamics of spinor condensates is driven be the competition of interaction energy and quadratic 
Zeeman effect. 

In our experiment, we prepare ^^Rb F = 2 atoms in the fully magnetized state in a transverse 
direction, such that the spin vector is (F) = 2ex in the co-rotating frame (see Fig.[2](a)). Neglecting 
the spatial degree of freedom, this state maximizes the mean-field energy and is therefore stationary 
at zero magnetic field, just as the anti-ferromagnetic ground state which minimizes the mean-field 
energy. However, analogous to an inverted pendulum, the fully magnetized state is stationary 
but energetically unstable, and thus an ideal starting point for dynamical instabilities and pattern 
formation. 

In order to study spatial pattern formation in the simplest non-trivial case, we have developed 
a reliable method to create extremely elongated Bose-Einstein condensates in the optical potential 
of a single focused laser beam, with an axial trap frequency of only cj^ = 27r x 0.8 Hz. The con- 
densates are approximately 800 /xm long and the calculated Thomas-Fermi radii perpendicular to 
the long axis are 2.9 /xm horizontally and 4.5 /xm vertically. In the transverse direction, the spin 

healing length = 'sj 2mgi{n) ^3.4 /xm is of the order of the size of the condensate, effectively 

suppressing transverse spin structure and creating a quasi- ID geometry. A sufficiently homoge- 

^Spin-2 condensates have an additional contribution ^2^|*5'oP where Sq = C+2C-2 — C+iC-i + Co/^- This contri- 
bution is included in our numerical calculations in the low-field case, however its influence is minor and we neglect it 
in this Letter for the sake of simplicity. 
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neous magnetic field, whose axial gradient is canceled with a remaining curvature of 60 mG cm"^, 
is aligned with the axis of the trap. Stern-Gerlach separation during a short time of flight allows 
us to image the individual spin densities n|CmP with a spatial resolution of about 3 fim. Image 
distortions caused by imperfect Stern-Gerlach fields are reversed by careful post-processing. 

We observe the spontaneous emergence of wave-like spin textures along the axis across a wide 
range of magnetic fields. Fig. [T] provides an impression of these patterns as they emerge, Fig. [2] 
shows a possible interpretation in terms of local spin orientation. One can clearly identify spin 
textures with well defined symmetry and wavelength for both low and high magnetic fields, while 
the structure is more complicated in an intermediate regime. At low magnetic field, where the 
interaction energy dominates (q <C gin), the mode is anti- symmetric with respect to positive and 
negative mp, indicating an alternating non-zero axial magnetization. Our data seems to suggest 
that both the wavelength and growth rate scale with the magnetic field. At high magnetic fields 
where the quadratic Zeeman energy far outweighs the interaction energy scale (q gin), mainly 
the = 0, ±1 spin components are modulated and the characteristic mode is symmetric with 
respect to = +1 and mi? = — 1. The axial component of the spin vector thus remains zero 
everywhere. The wavelength as well as the timescale of emergence of these patterns in the Zeeman 
regime do not change significantly as the quadratic Zeeman energy changes by a factor of four. 

Calculating cross-correlation functions of the individual spin densities allows to evaluate quan- 
titatively the periodicity of the modes as well as the dynamics of their growth (Fig. [3]). At 5 = 
0.25 G and B = 1.1 G, representative of the regimes q <C gin and q :» gin respectively, the 
wavelength is of the order of 27r x the spin healing length = 3A /xm in our experiment, namely 
28 /xm and 22 /xm respectively. The growth rate of the mode amplitude is half the measured growth 
rate of the correlation function, i.e. 20 s"^ and 28 respectively. In both cases the patterns sat- 
urate after about 200 ms. In contrast to classic pattern-forming system, the patterns developing in 
elongated spinor condensates are transient and decay into chaotic small-scale structure on the time 
scale of several hundred milliseconds. Whether this is an intrinsic effect driven by the same non- 
linear coupling that leads to the emergence of regular structure in the first place is an unresolved 
question. Other reasons for the observed decay may be the limited lifetime of atoms in the trap, 
leading to a loss of density, or additional mechanisms of interaction, such a dipolar forces 

We compare our observations to the Bogoliubov excitation spectrum of a homogeneous, one- 
dimensional spinor condensate prepared near the transversely magnetized state (Fig.|4]). It turns out 
that both in the interaction- and Zeeman-dominated regime singular unstable modes exist, indicated 
by an imaginary frequency over a certain range of wavevectors. The characteristic spin patterns of 
these modes indeed coincide very well with the patterns observed experimentally (Fig. [5]). Quan- 
titatively, the calculated wavelengths are 20% shorter than observed and the growth rates 40-65% 
higher. In light of the heavy approximations made in the calculation and the inherently limited 
validity of the Bogoliubov linearization, we consider this a reasonable agreement. 

^Note that the sHght residual axial magnetic field variation in our experiment could possibly influence the patterns 
observed at the lowest magnetic field B = SO mG. 
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Fig.[6]again illustrates the maximally unstable Bogoliubov modes by plotting their overlap with 
the Bloch sphere of all stretched states In the interaction-dominated regime, the mode can be 
qualitatively characterized (as already sketched in Fig.[2](b)) by a spin vector (F) rotating around 
the y-axis, with a full rotation corresponding to one spatial period of the mode. Thus, the axial 
and transverse magnetization oscillate out of phase to each other, on top of the fixed transverse 
magnetization of the initially prepared state. In the Zeeman-dominated regime, the evolution of 
the spin along the condensate axis is more complex and the local spin state generally lacks the 
rotational symmetry implied in Fig.|2](c). However it is still characterized by a plane of symmetry 
whose orientation rotates as the wave propagates. 

While linear stability analysis provides a satisfactory explanation of the variety and dynam- 
ics of modes observed in experiment, it is a purely classical approach in character. But our work 
also is a first step towards fundamental quantum mechanical aspects of ultracold matter. For ex- 
ample, having observed exponential growth, the question immediately arises what is the origin of 
the fluctuations triggering this growth A back-of-the-envelope calculation assuming 100% 
modulation after 200 ms and a growth rate of 20 s"^ provides an estimated modulation depth of 
2% initially. At our linear density of about 1.7 x 10^ cm"^, this corresponds to four times the shot- 
noise of the number of atoms contained within one wavelength of 28 /xm. Another interesting issue 
arises from the similarities between coherent dynamics in spinor BEC and optical non-linearities. 
Interpreting density patterns as interference between left- and right-propagating modes suggests 
the possibility of entanglement as in spontaneous parametric down-conversion. Finally, our results 
also demonstrate that purely local anti-ferromagnetic interactions suffice to generate ordered spin 
textures, in contrast to recent observations in ferromagnetic ^^Rb F = 1, where these textures were 
attributed to dipolar interactions 

Methods 

Experimental sequence Pre-cooled ^^Rb atoms in the F = 1, = —1 state are loaded from a 
magnetic trap into the single-beam dipole trap and evaporatively cooled to quantum degeneracy by 
ramping down the optical power over 13 s. During the dipole trap ramp, the magnetic field is set 
to its final value. The atoms are then prepared in the F = 2, mi? = —2 state by a microwave pulse 
and subsequently rotated into the fully transversely magnetized state by a radio-frequency pulse. 
After a variable time of evolution, the trap is switched off, the magnetic field is non-adiabatically 
switched to a high value and the Stern-Gerlach gradient field is applied. The total time of flight 
before absorption imaging is 12 ms. 

Image post-processing Linear and non-linear transformations are applied to the absorption im- 
ages to remove distortions introduced by the Stern-Gerlach procedure. The parameters of the 
transformations are independently determined from images obtained in the same way, but with 
regular patterns induced by a defined magnetic field gradient. Axial density profiles of individual 
rriF components are obtained by transverse integration of the absorption images. Correlation func- 
tions are then calculated in two different ways. The color-coded plots in Fig. [3] show normalized 
correlation functions of density profiles, i.e. at t = the overall density profile of the condensate 
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can be seen, and at Az = auto-correlations equal unity by definition. This method has proven 
useful for visually bringing out the emerging patterns. Growth rates in contrast are extracted from 
relative spin populations, in order to avoid the influence of the overall density profile. The cor- 
relation functions are not normalized to avoid artefacts. After spatial bandpass filtering centered 
around the observed periodicity, an rms amplitude is extracted and serves as a measure of the 
modulation contrast of the pattern. 

Calibration of parameters Transverse residual magnetic fields are carefully compensated to a 
level of 0.5 mG or less, while the axial magnetic field is accurately measured using rf spec- 
troscopy. The axial gradient of is zeroed at the center of the atomic cloud, with a remaining 
variation AB^ ~ 35 juG due to curvature at the ends of the cloud. Over the central 200 fim region 
displayed in Fig. [TJ the field-induced phase variation is less than tt for more than 200 ms. The 
average density of condensate atoms in the central region (n) = 6.5 x 10^^ cm"^ is obtained by 
comparing the global spin dynamics before the onset of pattern formation with the results of 

Linear stability analysis For low magnetic fields, we analytically linearize the homogeneous, 
one-dimensional, F = 2 Gross-Pitaevski equations using a Bogoliubov ansatz and assuming a 
stationary initial state. For a given set of parameters (gi,g2, q), a stationary state close to the fully 
transversely magnetized state is obtained numerically from the full non-linear equations of motion, 
and the corresponding Bogoliubov matrices are numerically diagonalized to obtain the spectrum. 
For high magnetic fields, we eliminate the fast dynamics induced by the quadratic Zeeman ef- 
fect by applying a unitary transformation U = exp{—iqF^t/h) and dropping all explicitly time 
dependend terms from the linearized equations. This approximation fully removes the magnetic 
field dependence and allows us to analytically find another stationary state close to the the fully 
transversely magnetized state. 
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Figure 1 Initial state (a) and saturated spin patterns (see text) arising in spinor conden- 
sates for increasing magnetic field (b-f). The individual 3D rendered false color absorption 
profiles in each image visualize the spatial density distribution for each of the five mag- 
netic projections = -2 ... + 2 of the F = 2 hyperfine manifold. Qualitatively different 
wave-like patterns arise both in the regimes of low and high magnetic field, which are 
separated by an irregular intermediate regime (d). Imperfections in the preparation and 
detection artefacts (e.g. optical interference fringes), which appear as fluctuations visible 
in the initial state (a), are distinctly different from the patterns of (b-f). The cases marked 
are representative of the initial state, the interaction dominated regime and the Zeeman 
dominated regime, respectively, and are referred to in the following figures. 
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Figure 2 Simplified illustration of the local spin vector in the Fig. [T] (a), (c) and (e) 
(marked '*'). (a) Initial state, the spin vector is fully stretched and points in the trans- 
verse X direction, (b) Interaction dominated regime, the spin vector is fully stretched but 
rotates from axial to transverse and back over one spatial period, (c) Zeeman dominated 
regime, the local spin vector has zero length but still a defined orientation [rriF = state) 
that rotates from axial to transverse and back. Compare also to the more realistic spin 
representation in Fig. [6} 
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Figure 3 Cross-correlations of spin densities, (a) Interaction dominated regime 0.25 G, 
(b) Zeeman dominated regime 1.1 G. Color-coded plots (left and right) of the normalized 
correlation as a function of relative axial displacement and time of evolution allow to de- 
termine the symmetry and spatial periodicity of the evolving modes. In (a), the mode is 
antisymmetric with respect to exchange of tuf = ±1 whereas the = component is 
not involved in the pattern formation. In (b), the mode is symmetric with respect to ± 1 
and antisymmetric with respect to = and both = ±1. The correlation contrast 
(center) can be tracked over two orders of magnitude and allows us to extract the growth 
rate of the mode and the identification of linear, saturation and decay regimes. 
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Figure 4 Bogoliubov spectra resulting from a linear stability analysis of a homogeneous 
system, (a) Interaction dominated regime 0.25 G, (b) Zeeman dominated regime ocG. For 
F = 2, five pairs of complex conjugate branches exist (of which only the positive real 
frequencies are plotted). In both cases the highest real frequency modes correspond to 
density waves (no spin modulation). Branches with positive imaginary frequency indicate 
unstalDle modes. The experimentally observed modes (crosses) are qualitatively repro- 
duced, but occur at slightly lower wave vector and have a lower growth rate than the most 
unstable modes (squares), which we attributed to the approximations in the model. 




Figure 5 Comparison of simulated Stern-Gerlach images with experimental data in the 
linear growth regime, (a) interaction dominated regime 0.25 G, (b) Zeeman dominated 
regime ocG. This figure demonstrates that the mode pattern of the calculated most un- 
stable mode actually conincides with with the experimentally observed ones during linear 
growth, despite quantitative deviations in wavelength and growth rate. 
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Figure 6 Bloch sphere projections of the calculated unstable mode, scaled to span the 
full color range from blue (minimum) to red (maximum), (a) Interaction dominated regime 
0.25 G, (b) Zeeman dominated regime oc G. For comparison, the Bloch sphere projections 
corresponding to Fig.|2]are also displayed (boxes on the left-hand side). These projections 
confirm the general features of the spin vector illustration in Fig. [2] but also show some 
differences. In the interaction dominated regime, compared to a fully stretched state 
rotating around the x axis, the actual mode rotates on a cone rather in the yz plane. The 
fully axially magnetized state is thus not reached in the linear regime. In the Zeeman 
dominated regime, the mode does not have the full ttif = Q symmetry and shows some 
additional modulations, however it is still characterized by a similar plane of symmetry. 
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